B — > XcO^K decays: a model estimation 
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> 1 Abstract 
oo ' 

CN ■ In this paper, we investigate the vertex corrections and spectator hard scattering contributions 

to B — > XcQ,2K decays, which has no leading contribution from naive factorization scheme. A 



non-zero binding energy b = 2m c — M is introduced to regularize the infrared divergence of the 
vertex part. The spectator diagrams also contain logarithmic and linear infrared divergences, for 
which we adopt a model dependent parametrization. If we neglect possible strong phases in the 
hard spectator contributions , we obtain a too small branching ratio for Xd)K while too large 
one for Xc2K, as can be seen from the ratio of the branching ratio of B + — > Xc2^ + to that of 
B + — ► XcqK + i which is predicted to be 2.15^!^ in our model, while experimentally it should 
be about 0.1 or even smaller. But a closer examination shows that, assuming large strong phases 
difference between the twist- 2 and twist-3 spectator terms, together with a slightly larger spectator 
infrared cutoff parameter A^, it is possible to accommodate the experimental data. This shows 
that, for B — > XcOfiK decays with no factorizable contributions, QCDF seems capable of producing 
decay rates close to experiments, in contrast to the B — > J/ipK decay which is dominated by the 
factorizable contributions. 
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I. INTRODUCTION 



Hadronic B decays attract a lot of attention because of its role in determining the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, extracting CP-violating angles and 
even revealing physics beyond the Standard Model (SM). However in most deep 
understanding on the strong dynamics in hadronic B decays is prerequisite for the above 
purposes. 

Phenomenologically the naive factorization ansatz (NF) jl], supported by color trans- 
parency argument Q], is widely used in hadronic two-body B decays. However the unphys- 
ical dependence of the decay amplitude on renormalization scale indicates a prominent role 
of QCD corrections to NF. In this respect, B — > Xcs^K decays are of special interest as 
these channels vanish in the approximation of NF, due to the spin-parity and vector current 
conservation. Therefore they provide a good opportunity to study the QCD corrections to 
NF. It was generally believed that the branching ratios of these channels should be quite 
small as the QCD corrections are either suppressed by strong coupling a s or AgcWmb- But 
BaBar [3] and Belle jj] have found a surprisingly large branching ratio of B + — ► Xd)K + 
decay, 

f (6.0 + ?i±l.l) x 1(T 4 (Belle) , 
B(B + -> XcoK + ) = 1 - L8 ^ 1 ' ' (1) 

I (2.7 ±0.7) x 1(T 4 (BaBar). 

Actually this large branching ratio is even comparable, for example, to that of B — > x c \K 

decay which is not forbidden in NF. Another surprising observation is that, the upper limit 

of B — > Xc2K decay is roughly an order of magnitude smaller than the observed branching 

ratio of B + — > XcoK + decay p, 

B(B + -> Xc2K + ) < 3.0 x 10" 5 (BaBar) , (2) 

while naively the branching ratios of B — > Xco,2K decays are expected to be at the same 
order. 

In the following we shall discuss these decay channels using the QCD factorization 
(QCDF) approach |6j. In this framework, the final state light meson is described by the 
light-cone distribution amplitude (LCDA), while for the P-wave charmonium Xc0 2, we shall 
adopt the covariant projection method of non-relativistic QCD p|. It is well known that, for 
the inclusive decay and production of P-wave charmonia, the color-octet mechanism must 
be introduced to guarantee the infrared safety. However it is still unclear how to incorporate 



this mechanism in a model-independent way into exclusive processes. Thus the decay am- 
plitudes A(B — > Xco,2K) would be inevitably infrared divergent when only the color-singlet 
picture is adopted for Xc, which is shown explicitly in $]. Thus strictly speaking, the QCDF 
approach is not applicable for B — > XcO^K decays due to the breakdown of factorization. 

In this paper, to get a model estimation, we will introduce the binding energy b = 2m c —M 
|t| as an effective cut-off to regularize the infrared divergence appearing in the diagrams of 
vertex corrections (see Fig. 1). In fact, the logarithmic divergence ln(fe) term in the limit 
6^0 for the vertex corrections in B — > Xca^K decays is similar to the ln(6) term found in 
Ref. p] for the production of P-wave charmonium in e + e~ collisions. As for the spectator 
scattering contributions, there appears logarithmic divergence at twist-2 level and linear 
divergence at twist-3 level. Phenomenologically we shall parameterize these divergence as 
lnfm^/A/j] and jtib/A/, respectively, where the non perturbative parameter = 500 MeV 
again acts as an effective cut-off to regularize the endpoint divergence According to 
the QCDF approach, all other contributions are power suppressed by AQCD/mb- 

We find that, with the above method, the branching ratio of B + — ► XcoK + decay is about 
0.78 x 10 -4 , which is several times smaller than the experimental measurements. At the same 
time, we also get the branching ratio of B + — > Xc2^ + decay at about 1.68 x 10~ 4 , which 
is significantly larger than the upper limit 3 x 10~ 5 observed by Babar [5j. But the above 
estimation is very crude in that the strong phases effects are completely ignored. Notice 
further that for the spectator contributions, there contains only logarithmic divergence at 
twist-2 level, while linear divergence appears at the twist-3 level, the strong phases of the 
twist-2 and twist-3 spectator terms could be quite different. We then briefly discuss the 
potential strong phases effects and argue that very different strong phases between twist-2 
and twist-3 spectator terms together with a slightly larger seems to be able to reproduce 
the experimental hierarchy B(B + — > XcoK + ) ^> B(B + — > Xc2K + ). 

II. VERTEX AND SPECTATOR CORRECTIONS 

In the QCDF approach, K meson is described by the following light-cone projection 
operator in momentum space 0] 
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FIG. 1: Order of a s contributions to B — > XcjK decay. Fig. (a)-(d) and Fig. (e)-(f) are called 
vertex corrections and spectator scattering diagrams, respectively. 



where I is the momentum of K meson and Zi^) is the momentum of quark (antiquark) in K 
meson. and $p(x) are leading twist and twist-3 distribution amplitudes of K meson, 

respectively. It is understood that only after the factor Z 2 • h i n the denominator is canceled, 
may we take the collinear limit l\ = xl, I2 = (1 — x)l. Notice that in principle we could also 



start directly from the original light-cone projector of K meson in coordinate space [ll|, and 
the physical results should be the same. But in this case care must be taken that, only with 
a proper regularization, can one do the relevant convolution integrals correctly. The readers 
may refer to the appendix of for further details. 

Since P-wave charmonium Xcj is involved, we shall use covariant projection method 
to calculate the decay amplitude 



q=0 



(4) 



Here A is the standard QCD amplitude for cc production, amputated of the heavy quark 
spinors, C\ = Sij/\/3 is the color singlet projector. While II" is the S = 1 heavy quark 
spinor projector 

^i=^^{^-i-Ai a i^ + i + A > ( 5 ) 



where P is the momentum of charmonium and 2q is the relative momentum between the cc 
pair in Xcj- ^aB * s polarization tensor of Xco,2 which satisfies the following sum over 
polarization relation 







c(o)c(o) _ - rr n 



1 



,(2)-(2) 1 — - „ „ , I 



^ad^a'P' ~ o (naa'Il/3/3/ + n Q , / g/il j a Q ,/) — -Xla/jITo/^/ , (6) 



with 

n_„ = o -+- : 

M 2 



n a/ ? - -# Q/ 3 + " / • (7) 



Here M is the mass of Xcj- 

For charmonium B decays, we shall start with the effective Hamiltonian ^| 

H ef f = ^ IvaVZiCtWM + C 2 (v)Q c M) - WEWw} 



(8) 



where Cj are Wilson coefficients which are perturbatively calculable and Qi,2(Q3-6) are the 
effective tree (QCD penguin) operators. Notice that we have dropped the electroweak pen- 
guin contributions here which are numerically negligible. The four-quark effective operators 
are defined as 

Ql — (<laba)v-A(C[3Cp) V -A , Q°2 = (s a bf 3 ) v _ A (c l3 C a ) V ^A , 

(9) 

Q3.5 = (s a b a )v-A XX<?/3<?/3)v=F^ j Q4.6 = (s^a)y-iE(?«9(3)yT^ • 

1 q 

Here q denotes all the active quarks at the scale \i = 0(rrib), i.e. q=u, d, s, c, b. While a 
and j3 are color indices. 

It is then straightforward to get the decay amplitude of B — > Xco^K decay by considering 
the vertex and spectator corrections drawn in Fig. 1, 

MB ^) = ^™ %%y<Mh - V a V i(Ci + O) x 

F o~* K ^/(o,2) + ^ /bS- (/(c3) + /(5S))^ ' ( 10 ) 

where i?i(0) is the derivative of the Xcj wave function at the origin and Fq^ k the form 
factor of B — > K. The function f 1 represents the contributions from vertex corrections while 
f II2 {f 113 ) arising from the twist- 2 (twist-3) spectator contributions. The vertex function f 1 
is actually infrared divergent and therefore depends on the binding energy b = 2m c — M. 
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In the following we shall keep ln(o/M) term and drop the terms suppressed by b/M. The 
explicit expressions of fL 2 ) are as follows 

~-b 



j _ 2m B ((l + 12a)(l-4a) + 16alii[4a]) 
~ (1 -4a) 2 VM 11 



M 



32£^*p a B p B ^((l + 12a)(l - 4a) + 16a In [4a]) 
f 2 = — m 



+" /Ofin + 0(b/M) , 
-b 



M 



+ fLn + 0(b/M) , (11) 



ms(l — 4a) 3 

where a = ml/ml , and /g n is the finite part of the function f 1 in the limit b/M — > 0. The 
explicit expressions of /g n are as follows, 



/ofin ~~ 



-m B 



| -6 - 221n2 + 4a(26 + (15 - 561n2)ln2 + 



2(l-4a) 2 (l-2a) 3 V3a 
8a 2 (65 + 52 In 2 - 84 In 2 2) + 384a 4 (l + 2 In 2) + 2a(-85 + 28 In 2(-l + 6 In 2)) + 

32a 3 (-23 - 32 In 2 + 14 In 2 2)) - 8 In a + 4a(-(l - 4a) 2 (5 - 24(1 - a)a) x 



In 



■1 + 4a 



+ 91na + 2(a(-3 + 4a)(13-46a + 56a 2 ) - 



4(1 - 2a) 3 In [64a]) In a + 16(1 - 2a) 3 In 2 In [-1 + 4a]) - 



A fin — 



64a(l - 2a) 3 ( Li2 
c(2)* n aJ> 

C a/3 PbPb 



2 -4a 



1 -4a 
32a 



+ Li2[l -4a] -Li2 



1 -2a 



)}■ 



Am By /a [(l-2a) 3 (l-4a) 

-8a(3a - l)(lna - In [4a - 1]) 



1 -4a. 

4 In 2(1 - 2a) 2 + (1 - 4a)(4a 2 (l + 2 In 2) - 1) 

V A B[a]\ , (12) 



where the function Vab[g] denotes the finite part of vertex corrections from Fig. (la-lb). 
The analytical form of Vab[o] is too complicated to be shown here, but numerically it has a 
very mild dependence on the parameter a, for example, 



V AB [0.1] = 11.3 , 



As for the spectator functions, we have 



V AB [0.15] = 11.9 



f ° = mt (l-4a)VS5 Jo ( ~T Jo dy — { - 8 " + (1 " 4a)v) ■ 

Jo 4 Jo y z 



(4a+(l-4a)y) , 



/o 



773 



m 6 (l - 4a) 2 V3a m fe 



J2 — 



m 3 (l — 4a) 3 Jo 

'(2)* 



dy- 



m%(l — 4a) 4 m& Jo £ Jo y 2 



Here £ is the momentum fraction of the light spectator quark in the B meson, and y = 1 — y 
the light-cone momentum fraction of the quark in the K meson which is from the spectator 
quark of B meson. Notice that our expressions for twist-2 spectator function /j^l) are 



consistent with those of 



III. NUMERICAL RESULTS AND DISCUSSION 



To get a numerical estimation on the branching ratios of B — > Xco^K decays, several pa- 
rameters appearing in Eqs. ()10I~I13|) should be first decided on. The derivative of Xcj wave 
'unction at the origin |-R'(0)| may be either estimated by QCD-motivated potential models 
1^ | . or extracted from Xcj decays 0. |-R'(0)| 2 varies from 0.075 GeV 5 to 0.131 GeV 5 in 
different potential models |lj] while using Xcj decays, for instance [l^ . 



x 36 4 2 \R'(0)\ 2 ( 16a, . 
T(Xc2 - 77) = t^L^f- ( o~~) . ( 14 ) 



c 



it is easy to get \R'{0)\ 2 = (0.062 ± 0.007) GeV 5 if we take m c = 1.5 GeV. This result is 
a little bit lower than, but still consistent with the potential model calculations, especially 
considering that it is very sensitive to the choice of charm quark mass. In this paper, we 
shall take ^'(O)] 2 = (0.10 ± 0.03) GeV 5 as input. 

For the binding energy, if we take m c = 1.5 GeV, the ratio b/M is about — 0. 1 1 ( — 0.16) 
for Xco(Xc2), while a = m 2 /ml ~ 0.1. The QCD scale \i should be order of ^m b A, as in 
charmless B decays, which is about (1 ~ 1.5) GeV. In the following we shall fix the scale 
/i = 1.3 GeV with a s = 0.36. Notice also that the Wilson coefficients should be evaluated 
at leading order, to be consistent with the leading order formula of Eq. (fTU|). 

d = 1.26 , C 4 = -0.049 , C 6 = -0.074 . (15) 

The relevant CKM parameters are chosen to be A = 0.83 and A = 0.224. 

As for the spectator contributions, we adopt the following LCDAs for the final K meson, 

My) = %(i - + E ^ /2 \2y - 1) , My) = i , (ie) 

n>l 



where C^^ 2 \x) are Gegenbauer polynomials. The parameters a n are set to be 17]: 

d\ = 0.17 , ct2 = 0.115 , a>4 = 0.015 , a 3 = a n> 4 = . (17) 
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Then logarithmic and linear divergences appear in Eq.([13p. which may be phenomeno- 
logically parameterized as 



dy m B 
— = In-— 

Aft 



dy 

,2 



m B 
Aft 



(18) 



y J y 

with Ah = 500 MeV. Notice that the above parametrization of linear divergence would 
violate the power counting of QCDF, but we do not have better way yet to deal with it. 
This is clearly a very rough estimation, for example, we do not consider here the strong 
phase effect. We also know little about B wave function, but fortunately only the following 
integral is involved which may be parameterized as 

MO 



1 



TUB 



(19) 



and we shall simply fix Ag = 350 MeV in our calculation. The chirally enhanced ratio tk = 
fix/rrib is chosen to be 0.43lo;og, which corresponds to taking m s (2 GeV) = (90 ± 20) MeV 



and (m u + (2 GeV) = 9 MeV. The form factor F ( 



B->K, 



in 



Xc 



may be read from |l7(, in 



which as stated, the uncertainty of form factor at q 2 ^ is likely to be smaller than that of 
q 2 = 0, which is about 12%. Therefore we will cite Fq~* k (rrv^J = 0.48 ± 0.06 as our input. 
The decay constants are set as fx = 160 MeV and fs = (210 ± 25) MeV. With the above 
input, we get 



B(B+ - X coK + ) 
B(B+ - Xc2 K + ) 



(0.78ig$) x 10- ZJ 



-0.69) >< 



'1 68 +0 ' 78N| x 1 ,v 1 



(20) 



We also show separately the contributions from vertex corrections and hard spectator scat- 
tering diagrams in Table I, with all the input parameters taken at their central values. For 



TABLE I: The numerical estimations of vertex corrections and hard spectator scattering contribu- 
tions, with all the parameters taken at their central values. The constant C = %h J§->k ■ 

c o 



Decay Channels 


f 


C * f 112 


C * f 113 


XcoK 


46.3 - 33.6i 


-43.1 


80.7 


Xc2K 


1.7+ 14.lt 


69.3 


68.3 



the case of XcqK channel, our results are approximately four times smaller than the average 
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of BaBar and Belle measurements, (3.0 ± 0.7) x 10~ 4 , while our prediction on B — > Xc^K 
decay is obviously too large compared with the experimental upper limit, 3.0 x 10 -5 . This is 
a little bit surprising, because for charmonium B decays, the theoretical results are normally 
a few times smaller than the experimental measurements. 

The careful reader may have noticed that in the above analysis we did not consider the 
uncertainty related to the parameter a = m^jvn^. In fact a larger a could enhance the 
branching ratio of B — > XcoK decay significantly, but unfortunately it would also enhance 
that of B — > Xc2K decay with similar magnitude. Notice that B — > XcapK share many 
common inputs, the ratio of branching ratios of these two channels should have mild depen- 
dence on the input parameters, for example it is independent on the parameter |i?'(0)|. Our 
numerical analysis shows that this is indeed the case, with a = 0.10±0.03: 

-1 c+0.26 +0.33 +0.36 +0.30 /91 \ 

K ~ B (B+ - Xc0 K+) ~ Z --°- 55 -°- 31 -°- 28 "°- 31 ' [Zl) 

where the uncertainties arise from the parameters a, rx, Fq~ >k and respectively. The 
above ratio is clearly in strong contradiction with the experimental hierarchy 72. 0.1 <C 1. 

Notice that the chirally enhanced power corrections, namely twist-3 spectator contribu- 
tions in this case, have been included in the above estimation. For the rest part of power 
corrections, there is no systematic way to estimate them yet. But since the power correc- 
tions are suppressed by Aqcd/ 777 *; intuitively they might lead to an uncertainty of about 
20% to the decay amplitude, which is unlikely to be able to change our estimation Eq. (f2~Tj) 
dramatically. 

In our model, the parameters and Xb will introduce additional uncertainties to B — > 
Xco,2K decays. It is very unlikely that we could reproduce the experimental observations 
by fine tuning A#, because although a larger A# would lead to a smaller branching ratio 
for Xc2K decay, it would also make the already too small branching ratio of XcqK decay 
even smaller. However a larger does help to close the gap between our predictions and 
the experimental data, due to the fact that a larger A/j will lead to a significantly smaller 
branching ratio for Xc2^ decay while XcoK decay does not change much. Of course we can 
not choose a too large A^, say larger than 1 GeV, because it is anyway a non perturbative 
parameter. As an illustration, we take A^ = 700 MeV and get 

B(B + -> XcoK + ) = 0.78 x 10~ 4 , B(B + -> Xc2K + ) = 0.74 x 10" 4 . (22) 
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Although it seems to be on the right way, this effort alone is still not enough to accommodate 
the experimental data. 

Let us take a closer look at the decay amplitudes. From Table I, it is clear that the spec- 
tator hard scattering mechanism is dominant in B — > Xc2^ decay and also very important 
f° r XcoK channel. Furthermore there is significant destructive (constructive) interference 
between the twist-2 spectator term and the twist-3 one for XcoK (x C 2K) mode. This is 
probably the reason that we get too small XcoK decay as well as too large Xc2K decay in 
our model. Notice that there are logarithmic and linear divergences appear in the spectator 
contributions, which are parameterized by Eq. ()18j) . It is obviously a very rough model es- 
timation and for example, strong phases effects are completely ignored. It is also reasonable 
to assume that the strong phase of twist-2 spectator term could be different from that of 
twist-3 part. As an illustration, the endpoint divergences could be parameterized as p: 

/^ =ln ^ (iw »,,), 1%^+^), (23) 

J y A h J y 2 Aft 

with < p < 1 and the phase 6 completely free. In the above equations, (p2,#2) denotes 
the parameters for twist-2 spectator term and (p3,# 3 ) for twist-3 one. In this case, the 
interference effects and therefore the predictions of the branching ratios, could be changed 
dramatically. For example, if we take a somewhat extreme case 

p 2 = 0.6, #2 = vr, p 3 = 0, #3 = 0, 

with A/j = 600 MeV while keep all other input parameters fixed at their central values, we 
will get 

B(B + - XcoK + ) = 3.3 x lO" 4 , B{B + - Xc2 ^ + ) = 1.7 x 10~ 5 , (24) 

which are in good agreement with the experimental observations. Certainly, due to the non 
perturbative nature of the above strong phases, there is strong model dependence of our 
predictions. Therefore it is not so meaningful to fine tune the parameters to get the best 
fit of the experimental data. The key point here is that, different strong phases between 
twist-2 and twist-3 spectator terms might be able to account for the experimental hierarchy 
that B{B — > XcoK) is at least an order of magnitude larger than B{B — > Xc2K). 

The authors of ref. also studied the B — > XcqK decay with the same QCDF method. 
But they used the gluon mass and gluon momentum cutoff, instead of binding energy adopted 
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in this paper, to regularize the infrared divergences of the vertex corrections. Another 
difference is that they calculated the spectator contributions directly from the original light- 
cone projector of K meson in coordinate space and got a different result from this paper. 
They claimed that the difference was due to the light-cone projector adopted in this paper 
which is inappropriate for XcjK channels: to get the projector Eq.fjSJ) from the original one in 



coordinate space the integration by parts has been used and the boundary terms were 
dropped. However because of the linear singularities appeared in the above calculations, 
the boundary terms seem to be divergent and thus the justification of using the integration 
by parts is in doubt in this case. But Beneke has elaborated on this subtle point in the 
appendix of ref. [l2j and it is shown there that the boundary terms are indeed zero provided 
the propagators are regularized carefully when they go close to the mass-shell. Therefore 
the integration by parts can be used here and the light-cone projector adopted in this paper 
is justified. Certainly, with a proper regularization, the calculation starting directly from 
the coordinate space projector should give the same results as this paper. 

Most recently, the B — ► XcoK decay was discussed by using the PQCD method |22| . 
Notice that the vertex corrections were not included in their calculations, and the spectator 
contributions alone are enough in their paper to account for the experimental data. It would 
be very interesting to see whether they could also reproduce the very small branching ratio 
of Xc2K channel observed by Babar, which has not been done yet. 

The B — > XcoK decay was also analyzed with light-cone sum rules Although 
there are some discrepancies in their papers, they agreed on the point that their results 
were too small to accommodate the experimental data. A large charmed meson rescattering 
effects B — > £)(*)£)(*) — > XdopK could account for the surprisingly large B — > XcoK decay 211 ]. 
but generally it will also lead to a large branching ratio for Xc2K mode. 

In summary, we discuss in this paper the vertex corrections and spectator hard scatter- 
ing contributions to B — > Xco^K decays. Since there is no model independent way yet to 
estimate the color-octet contribution to exclusive processes, it is no wonder that the vertex 
corrections here are infrared divergent. The non-zero binding energy b = 2m c — M XcJ makes 
the charm quark slightly off-shell inside Xcj, and effectively acts as a cut-off to regularize the 
vertex part. There are also less serious logarithmic and linear endpoint divergences which 
appears in the spectator contributions and are parameterized in a model-dependent way as 
usually done in charmless B decays. This means that the spectator diagrams are actually 
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dominated by soft gluon exchange, which in a sense could be viewed as a model estimation 
of color-octet contributions. Then our numerical analysis predicts the branching ratio of 
B + — > XcoK + decay to be about 0.78 x 10~ 4 , about four times smaller than the experimen- 
tal observations, while for B + — ► Xc2^ + decay, we get 1.68 x 10~ 4 , which is about five times 
larger than the experimental upper limit. But concerning the large theoretical uncertainties, 
it is more interesting to consider the ratio R = B(B + — > Xc2K + )/B(B + — > XcoK + ), in which 
a large part of the theoretical uncertainty can be eliminated. Numerically we find the ratio 
to be R = 2.15l{j;?j!, in sharp contrast to the experimental observation that this ratio should 
be about 0.1 or even smaller, if the BaBar analysis of the upper limit of B(B + — > Xc2K + ) 
decay will be confirmed by farther measurements. We then have a closer look at the decay 
amplitudes. One observation is that, XcqK channel is not very sensitive to the spectator 
infrared cutoff parameter while a larger could reduce the branching ratio of Xc2K 
decay significantly. Another observation is that, in our model there is large destructive 
(constructive) interference between the twist-2 and twist-3 spectator terms for XcoK (x&K) 
mode. But notice that the twist-2 spectator contributions contain only logarithmic endpoint 
divergence while twist-3 ones contain more severe linear endpoint divergence, it should be 
reasonable to assume that their strong phases could be quite different. Since the interfer- 
ence effects are very sensitive to the strong phases difference, this might change our model 
predictions dramatically. As an illustration, we then show in an explicit case that, with a 
slightly larger A^ and large strong phases difference between twist-2 and twist-3 spectator 
terms, our predictions are in good agreement with the experimental data. 

In conclusion, what we have shown in this paper, is that by adjusting the parameters 
for the spectator hard scattering contributions, as with the annihilation terms for charmless 
B decays, QCDF is able to produce appreciable non factorizable contributions to B + — ► 
Xco,2K + decays close to experiments, in contrast with the B + — > J / ipK + decay which needs 
a large factorizable contribution in addition to the small non factorizable one obtained in 
QCDF H- 
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